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Abstract 

Antineutrinos born inside the Earth ( "geoneutrinos" ) carry out in- 
formation of fundamental importance for understanding of the origin 
and evolution of our planet. We show that Baksan Neutrino Observa- 
tory is one of the best sites for detection and analysis of geoneutrinos 
using large liquid scintillation spectrometer. Also we present a short 
story of concept of Earth as antineutrino source (1960 - 2004 yy). 

Introduction 

In this paper we consider future studies at BNO (Baksan Neutrino Obser- 
vatory of Institute for Nuclear Research RAS) of terrestrial antineutrinos 
z/ e coming from beta decay of Uranium and Thorium daughter products 
( "Geoneutrinos" ) . Geoneutrinos bring information on the Uranium and Tho- 
rium abundances and radiogenic heat sources inside the Earth, which are of 
a key importance for understanding of the formation and subsequent evolu- 
tion of our planet [1] . Geoneutrino is a part of future studies at BNO of low 
energy z7 e fluxes of natural origin aimed for obtaining information on their 



sources, which is otherwise inaccessible. For this purpose a large target mass 
z/ e scintillation spectrometer is planned with the inverse beta decay 

v e + p — > n + e + (1) 

as the detection reaction. 

The program also includes: 

• Estimation of frequency of gravitational collapses in the Universe by 
detection of isotropic flux of z7 e , [2, 3]. 

• Test of the hypothesis that a self-sustaining nuclear chain reaction is 
burning at the center of the Earth ( "Georeactor" ) [4—7]. 

• Search for solar antineutrinos, which can be produced through 8 B neu- 
trino spin-flavor precession. 

• We note also that no spectroscopic information has so far been obtained 
on the solar pep, 7 Be and CNO neutrinos (i/ e ). We hope that in this 
important field progress can be achieved at BNO with u e , e scattering 
as detection reaction. 

• Quite recently old ideas [8] on using detection for nuclear reactor control 
(Plutonium production, nonproliferation etc.) have been revived: (See 
"Neutrino and Arms Control" 2004, http://www.phys. hawaii.edu~jgl/ 
nacw.html). BNO can join international efforts in this field. 

KamLAND Collaboration in Japan, using reaction (1) has recently demon- 
strated revolutionary progress in i7 e detection technique. KamLAND has 
already reported that a few geoneutrino events have been observed [9] and 
moves towards detection of soft solar neutrinos. BOREXINO Collaboration 
at LNGS, Italy, is planning to start soon experiments in the same direction. 
Under the circumstances question can be asked: Is it necessary to start new 
projects in the same field? 

As far as geoneutrinos are concerned, the answer is that their flux is 
expected to depend strongly on the point of observation. The intensity of 
z/ e is predicted to be maximal at high mountains sites (Himalaya, Caucasus 
where BNO is situated), to be minimal in the ocean far from continents 
(Hawaii, and may be Curacao [7]) and have intermediate values in Japan 
and Italy. 
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We note also that for success of the experiments listed above the level 
of the backgrounds is of a decisive importance. To a great degree the back- 
ground is produced by cosmic muons and by v e flux from power reactors. 
At BNO located at the 4800 mwe rock overburden, the muon flux is much 
lower than in KamLAND and BOREXINO experiments (Fig. 1). The flux 
of reactor u e at BNO is also 10 and 2 times lower than it is at mentioned 
laboratory sites respectively [6]. 

We conclude that BNO is one of the most promising sites to build a mas- 
sive antineutrino scintillation spectrometer for these studies. Essentially we 
propose to prepare the next step in low energy antineutrino and neutrino 
physics, which, using experience accumulated in KamLAND and BOREX- 
INO projects, will provide further progress in the fields discussed here. Im- 
portance of this goal is confirmed by recently published idea of LENA exper- 
iment [10]. 

In the next sections we give a short overview of the history of geoneutrino 
problem, which is now almost 45 years old, next we consider geoneutrino 
models and will finish with the results we hope to obtain at BNO. 

1 Geoneutrinos: First Estimations of Fluxes 
and Search for Detection Reaction 
(1960-1984) 

Marx and Menyard [11] in 1960 were the first to point out that the Earth is 
a source of antineutrinos, which are produced in U and Th daughter product 
beta decay, in decay of 40 K and of some other long-lived nuclei. They assumed 
that natural radioactivity is concentrated in a thin upper layer of the Earth 
and estimated the v e flux at the surface: 

F ~ 6.7 x 10 6 cm~V\ (2) 

which is not far from modern values. Radioactivity inside the Earth is un- 
known. If the concentration were not to decrease with the depth, the flux 
could be hundred times higher. 

"In the far future an experiment will be needed to set an upper limit 
on the antineutrino activity of the Earth. This probably is the only way to 
get information on the composition of the substance in the deep layers of the 
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Earth" wrote M. A. Markov in 1964 [12]. Markov, as far as we know, was the 
first to mention the process (1) as a possible geoneutrino detection reaction. 
"Because of high threshold, — he continued, — the number of active v e in 
this case is very low" . He stressed also the importance of finding a detection 
reaction with lower energy threshold. (Typical geoneutrino energy spectrum 
is shown in Fig. 2.) 

Two years later Eder [13] explained the gradual slowing down of the 
circulation rate of the Earth and the increase of its radius by a large amount 
of radiogenic heat sources in the mantel, "sufficient to blow up our planet" . 
The relevant geoneutrino flux was found two orders of magnitude higher than 
(2). Eder also mentioned process (1) as a possible detection reaction. No 
details of detector layout and detection principle were discussed at that early 
time. 

Antineutrinos of sufficiently low energy can, principally, be detected via 
resonant reactions of induced orbital electron capture in otherwise stable 
nucleus A(Z) [14]: 

u e + e- + A(Z) -> A(Z-l). (3) 

The reaction signature can be appearance of radioactive nucleus A(Z — 1) 
in the irradiated target (radiochemical method). Antineutrinos are captured 
only in a narrow energy band near the resonant energy E res m T max , where 
Tmax is maximal kinetic energy of the electrons in the decay of AiZ — 1): 
A(Z — 1) — > A(Z) + e~ + v e + T. Clearly, the resonant energy E res is 
2mc 2 = 1.02 MeV lower than the threshold energy in the "usual" inverse 
beta decay reaction on the same nucleus: 

v e + A(Z) -> A(Z - 1) + e+. (3.1) 

The transition probability w per unite time between initial and final states 
in reaction (3) can be written as: 

w ~ (ft)' 1 | V(0) | 2 dF(E res )/dE. (4) 

Here | ip(0) | 2 is the probability density for atomic electron at nuclear 
surface, ft is the reduced lifetime of beta decay (3.1) and dF(E res ) / dE 
cm -2 s" 1 MeV -1 is the spectral density of incoming v e flux at resonant energy. 
The density | ip(0) | 2 rapidly increases with atomic number, approximately 
as Z 3 , the spectral density is high at energies lower than 1.8 MeV (Fig. 2), 



4 



the inverse beta decay reaction (3.1) at high Z is suppressed by positron 
repulsion in the final state. As a result the resonant capture dominates over 
(3.1) for sufficiently heavy target nuclei. 

Potentially, the resonant capture reaction, unlike reaction (1), can search 
for terrestrial 40 K antineutrinos, which is of a vital importance for geophysics. 

In 1968 y Markov and Zatsepin organized in Moscow an International 
Workshop on Neutrino Physics and Astrophysics, which later transformed to 
regular "NEUTRINO" Conferences. At this Workshop R. Davis considered 
for geoneutrino detection purposes two resonant reactions: 

35 C1 ^35g ; Ereg = Q 1? Me y 209 Bi ^209p b) R ^ = q g 4 Me y ^ 

with subsequent radiochemical separation of radioactive 35 S and 209 Pb. 

Marx ([15], 1969) considered various models of terrestrial U, Th and K 
distribution, estimated relevant z/ e fluxes and recognized the importance of 
both resonant and inverse beta decay reactions for geoneutrino detection. 

Krauss, Glashow and Schramm ([16], 1984) calculated radiogenic heat 
powers, z/ e spectra and fluxes due to U, Th and K (Table l).They assumed 
spherically symmetric distribution of radioactive isotopes in a 100 km thick 
outer layer of the Earth; thus their results correspond to some average point 
on the surface. Most attention in [16] is given to geoneutrino registration via 
resonant and inverse beta decay reactions; the cross-sections for large number 
of nuclei have been calculated. In the most favorable cases the transition 
probability (4) was found to be: 

w ~ (3 - 6) x 10~ 3 7year. (6) 

This is several thousand times lower than is observed in solar Ga^ Ge ex- 
periments. . . "Detection of terrestrial antineutrinos will require sophisticated 
new technology. We can envision several generations of experiments aimed 
at measuring different parts of the antineutrino spectrum" . The first of these 
generations based on reaction (1) has surpassed in 2002 y the sensitivity level 
(6) and started detection of U and Th antineutrinos. Potassium antineutrinos 
still wait for their opportunity. 
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2 Radiogenic Heat - and Antineutrinos 
Models 

Here we widely use information taken from book [1], encyclopedias, from 
Hofmann papers [17] and Fiorentiny group [18] publications. Modern models 
distribute the Earth's antineutrino- and radiogenic heat sources U, Th and 
K between the crust and the mantle. The Earth's structure schematically 
shown in Fig. 3 includes components: 

• The crust. Thickness of the continental crust varies from ~70 km (at 
mountain sites like Hymalaya, Baksan etc.) to about 10 km, averaging 
around ~ 35 km. The continental crust contains a considerable part 
of the Earth's radioactivity. The oceanic crust is 5-6 km thick, covers 
about 3/4 of the Earth's surface; the U, Th and K concentrations are 
much lower than in the continental crust. 

• The outer frontier of the Earth's core lies at the depth of 2900 km. It 
is believed to be in a molten state because transverse seismic waves do 
not penetrate inside it. Inside the molten core is (partially) crystallized 
inner core. It is assumed that both cores are built up of iron, nickel 
and some amount of lighter elements; it is believed (by the majority of 
geochemists) that there is no radioactivity in the core, because U, Th 
and K are chemically "incompatible" with the main components. 

• The mantle, which is between the crust and molten core, is subdivided 
in two layers: the upper and the low. Volcanic outflows coming from the 
upper mantle show, that this layer is strongly "depleted" of radioactive 
elements. The lower mantle seems to be not so much depleted, but no 
systematic information is available on the subject. 

How much of U, Th and K is contained in the Earth (crust + mantle)? 
It is believed that meteorites of a special type (carbonaceous chondrites) 
originate from the same substance that formed the primordial Earth 4.55 
billion years ego. Starting from this point geochemists deduce the result: 

a(U) = 2 x 10- 8 g/g, a(Th)/a(U) = 3.9, (7) 

where a(U) and a(Th) are element mass abundances. 
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Using (7) one can find total mass of Uranium M(U) and total mass of 
Thorium M(Th), contained in the crust + mantle and can calculate relevant 
radiogenic heat powers H(XJ) and H(Th): 

Crust + Mantle : M(U) = 0.8 x 10 17 kg, M(Th) = 3.15 x 10 17 kg, 
Crust + Mantle : H(XJ) = 7.6 TW, H (Th) = 8.3 TW (8) 

Some features of U and Th decay are given in Tables 2, 3. The Earth in 
this model radiates L u = 1.1 x 10 25 v e j s emitted in U- and Th decay chains, 
which carry away 1.1 TW. 

The outlined model (sometimes called Bulk Silicate Earth, BSE) is shared 
by a great majority of geochemists and geophysicists. 

However, the BSE prescription does not work with K. Average K to U 
concentration ratio in carbonaceous chondrites is ~ 6 higher than is typical 
for terrestrial samples. This problem has been debated for decades. Two 
ideas have been proposed to explain K deficit in Earth's crust + mantle. 

1. Tiny particles of which the Earth was formed 4.55 billion years ago had 
at that time higher temperatures than chondrites and had lost most 
part of original K through evaporation. In this case K contribution to 
the radiogenic heat is H(K) ~ 2.5 TW. 

2. Actually there is no K deficit at all; just the major part of primordial K 
has gone into the core. In this case H(K) ~ 16 TW, most part of this 
power comes from the core. Recent laboratory experiments provide 
evidence that actually large quantity of K could enter the iron rich 
Earth's core [19]. 

Questioned is the theory that there is no U in the Earth's core. As already 
mentioned a self-sustaining fission reaction is supposed to burn there with a 
power release of 3—9 TW [4]. Also questionable is the type of meteorites, 
which ought to be used as a starting point for finding the amount of Earth's 
radioactivity [20]. 

In the BSE model the total radiogenic power H ra d = H(XJ) + H(Th) + 
H(K) ~ 19 TW is slightly lower than one half of the total Earth heat outflow 
(Htot ~ 40 — 42 TW). As an upper limit sometimes is considered a "fully 
radiogenic" model with U etc. masses increased by a factor ~ 2.2 so that 
H ra d — H tot . 
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3 Strategy of Research, Geoneutrino Fluxes, 
Observables 



Progress in the low energy v e detection technique (the CHOOZ experiment, 
KamLAND and BOREXINO proposals) opened practical opportunities to 
detect geoneutrinos. In 1997-1998 yy Rothschild, Chen and Calaprice [21] 
and Raghavan et al. [22] proposed to use large target mass liquid scintillation 
detectors to investigate the geoneutrinos with the inverse beta decay (1) as 
a detection reaction. 

Authors of [21] recognized that measurements are to be done at a number 
of geographical positions: 

• in mountain regions where contribution of core radioactivity is domi- 
nant and geoneutrino intensity is expected to be maximal, 

• in the ocean far from continents with dominant contribution from the 
mantle and minimal expected geoneutrion flux and 

• at Kamioka site in Japan and at LNGS in Italy with an intermediate 
level of geoneutrino intensity. 

Since 2002 y Fiorentini group in Italy has started extensive calculations 
of expected geoneutrino fluxes, detection rates and their dependence on the 
point of observation on the Earth's surface (see e.g. [18]). 

In scintillation detectors geoneutrino energy spectrum at E > 1.806 MeV 
is measured using delayed coincidences between spatially correlated positron 
and neutron produced in reaction (1). T ps E — 1.806 gives relation between 
positron kinetic energy T and the energy of incoming v e . In large target 
mass detectors the positron annihilation quanta are absorbed in the feducial 
volume and positron energy release E vis is : 

E vis m E — 0.8 (MeV), (9) 

Thus the minimal positron energy release at the reaction threshold when 
T = is E vis(mm) = 1.02 MeV. 

The numbers of Uranium iV(U) and Thorium iV(Th) geoneutrino events 
produced in the scintillator target exposed to incident z/ e fluxes F(XJ) and 
F(Th) are given by: 
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iV(U)/year 10 32 H = 13.2 F(U) x 10" 6 cm^s -1 

N(Th) /year 10 32 H = 4.1 F(Th) x 10~ 6 cm^s- 1 (10) 

(1 160 ton of CH2 based scintillator contains 10 32 H atoms) 

Geoneutrino fluxes F(U) from Uranium decay calculated in BSE model 
in [18] and [21] for a number of sites are shown in Fig. 4. In [18] higher 
values of Uranium abundances in the continental crust are assumed than 
those used in [21]. As a result differences in calculated fluxes for the same 
point of observation are rather large, which demonstrates good sensitivity of 
the neutrino method to the assumption on the distribution of radioactivity 
in the Earth. In Fig. 4 one can also clearly see a rapid decrease of the flux 
when the point of observation shifts from the mountains to the oceanic site. 

The total numbers of the Uranium and Thorium events iV(U) + iV(Th) 
(Table 4) also demonstrate large scattering of predictions found in [18,21,22]. 
The expected number of geoneutriono events at BNO near the Elbrus moun- 
tain is close to that found for geographical maximum in Himalayas. As can 
be seen in Fig. 5 Uranium and Thorium events can be separated and thus 
the ratio of U/Th abundances measured. One also can see that at BNO the 
background coming from nuclear power reactors is sufficiently low. 

The size and construction of the future detector and backgrounds require 
special consideration. Partially these issues have been considered in our 
recent publication on georeactor [6]. Good passive shield should protect 
the scintillator from the natural radioactivity of the surrounding rock, from 
PhM's glass and metallic support constructions. The anticoincidence system 
should effectively tag cosmic muons and showers. As already mentioned 
in the Introduction the muon flux at BNO is considerably lower than at 
Kamioka and LNGS (Fig. 1). Experience accumulated in KamLAND and 
BOREXINO projects shows that deepest purification of liquid scintillator 
can be achieved and extremely low level of radioactive contaminations can 
be reached, such as ~ 3 x 10~ 18 U g/g. 

In [6] we have considered LS volume sufficiently large (10 32 H atoms, 
1.16 kton) to confirm or reject theory that a nuclear reactor is burning at 
the center of the Earth. The geoneutrino studies seem to require larger tar- 
get mass detectors, may be as large as (1 — 2) x 10 33 H atoms. Indeed, (a) 
the detection rates due to oscillation effect will be ~ two times lower than 
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those given in the Table 4; (b) good statistics is needed to separate U- and 
Th geoneutrinos (Fig. 5) and, finally, (c) we consider here a new genera- 
tion experiment, which will operate after KamLAND. With an exposition 
(50 - 100) x 10 32 H • year it would be possible to accumulate at BNO 2000- 
4000 geoneutrino events. With statistics like this one can try to investigate 
the angular distribution of incoming geoneutrinos, which can give valuable 
information on the location of U and Th masses in the interior of the Earth. 

Conclusion 

Uranium and Thorium distribution in the Earh's interior and relevant ra- 
diogenic heat can be found via geoneutrino observations at a number of 
geographical points. 

The Baksan Neutrino Observatory is one of the most promising sites 
to build a massive antineutrino scintillation spectrometer for geoneutrino 
studies and also for investigations on low energy v e and v e fluxes of natural 
origin aimed for obtaining information on their sources, which are otherwise 
inaccessible. 

In 1984 y Krauss, Glashow, and Schramm wrote: '"Neutrino and antineu- 
trino astronomy and geophysics can open vast new windows for exploration 
above us and below". The first window is now opening at Kamioka. New 
windows can be opened at Baksan and elsewhere. Detector for Potassium 
antineutrinos is still to be invented. 
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Table 1: Potassium, Thorium and Uranium in the Earth's lithosphere: 
Masses (M), fluxes of u e (F) and radiogenic heat powers (if) according 
to [16] 



Element 


M 


F 


H 




10 17 kg 


10 6 cm^s" 1 


10 12 W 


U 


0.80 


3.5 


9.6 


Th 


3.8 


3.5 


7.3 


K 


5.2 x 10 3 


11 


1.8 




Total: 


18 


19 



Table 2: Energy balance (MeV/decay) in equilibrium 238 U and 232 Th decay 
chains: Total energy E tot , thermal (radiogenic) energy Eh and energy carried 
away by antineutrinos E v 





E t ot E u E v 


a38 U= ^ ^Pb + 8 4 He + 6z7 e 
232 Th ^ 2 08pb + 6 4 He + 4^ e 


51.7 47.7 4.0 
42.7 40.0 2.3 



Table 3: Specific radiogenic heat powers H, neutrino luminosities L v and 
radiated antineutrino powers W v 





#,TW/10 17 kg L u , 10 24 s-Vl0 17 kg W„,TW/10 17 kg 


238 TJ 
232 Th 


9.5 7.6 0.81 

2.6 1.6 0.15 



Table 4: Expected rates of U + Th events at various sites (exposition 
10 32 H • year, no oscillations) 



Site 


Rothschild et al* 


Raghavan et al. 


Fiorentini et al. 




[21] 


[22] 


[18] 


Himalaya 


65 




112 


Baksan 






91 


Gran Sasso 


53 


la: 134; lb: 50 


71 


Kamioka 


48 


la: 75; lb: 50 


61 


Hawaii 


27 




22 



* Calculated using Eqs. (10) and v e Fluxes from Ref. [21] 
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Figure 1. Muon flux underground (from V. Gavrin's talk at NANP-2003, 
http://www.nanp.ru/2003/program.html ) 



12 



U+Th+K 




12 3 4 



Antineutrino energy E, MeV 



Figure 2. Typical U+ Th + K geoneutrino energy spectrum. The arrow shows 
threshold of reaction (1) 




Figure 3. Internal structure of the Earth (schematic) 

1 - crust , 2 - upper mantle, 3 - lower mantle, 4 - outer molten core, 4 inner (solid) 
core. 




6 2 1 

Figure 4. Calculated Uranium geoneutrino flux F (10 cm" s" ) at different sites. 
1 - Fiorentini et al. [18], 2 - Rothschild et al. [21]. 




Figure 5 Expected geoneutrino-induced positron spectrum at Baksan, 

(exposition 10 32 Hyear, no oscillations) 

Solid line - U + Th, 1 - Th, 2 - signal from power reactors. 



